Structure maps for hep metals from first principles calculations 
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The ability to predict the existence and crystal type of ordered structures of materials from their 
components is a major challenge of current materials research. Empirical methods use experimental 
data to construct structure maps and make predictions based on clustering of simple physical param- 
eters. Their usefulness depends on the availability of reliable data over the entire parameter space. 
Recent development of high throughput methods opens the possibility to enhance these empirical 
structure maps by ab initio calculations in regions of the parameter space where the experimental 
evidence is lacking or not well characterized. In this paper we construct enhanced maps for the bi- 
nary alloys of hep metals, where the experimental data leaves large regions of poorly characterized 
systems believed to be phase-separating. In these enhanced maps, the clusters of non-compound 
forming systems are much smaller than indicated by the empirical results alone. 



INTRODUCTION 

Predicting the stable structures of alloys from their 
constituents is a major challenge of current materials 
research [l], E|. The traditional approach to this prob- 
lem is to extract trends from systems for which exper- 
imental data is available and apply them to uncharac- 
terized systems Q. These rules usually depend on sim- 
ple parameters, e.g. atomic number, atomic radius, elec- 
tronegativity, ionization energy, melting temperature or 
enthalpy, etc. Several well-known methods include the 
Hume-Rothery rules Q , the Miedema formation enthalpy 
method the Zunger pseudo-potential radii maps Q 
and Pettifor maps (zTQ- 

These methods produce results largely consistent with 
existing experimental data and have helped direct a few 
successful searches for unobserved compounds [9]. How- 
ever, they offer a limited response to the challenge of 
identifying new compounds because they rely on the ex- 
istence of consistent and reliable information for systems 
spanning most of the relevant parameter space. In many 
cases, reliable information is missing in a large portion 
of the parameter space — less than 50% of the binary sys- 
tems have been satisfactorily characterized This 
leaves considerable gaps in the empirical structure maps 
and reduces their predictive ability. 

High-throughput (HT) calculations of material prop- 
erties based on density functional theory (DFT) have 
been developed for theoretically guided material discov- 
ery and design In this approach the phase sta- 
bility landscape of alloys is explored by calculating the 
formation enthalpies of a large number of possible struc- 
tures, identifying the minima at various component con- 



centrations. A minimum-free-energy convex hull is con- 
structed from these minimum energy structures. These 
calculations give insights into trends in alloy properties 
and indicate possible existence of hitherto unobserved 
compounds. Recent advancements in this method, by 
efficiently covering extensive lists of candidate structure 
types [18l421|. make it possible to complement sparse ex- 
perimental data with ab initio total energy assessments. 
This development was envisioned by Pettifor a few years 
ago and here we present its first realization for binary 
metallic systems. 

We apply the HT approach to a comprehensive screen- 
ing of the 105 intermetallic binary systems of the hep 
metals. These systems have been generally less exten- 
sively studied than the cubic (fee and bec) metals, both 
experimentally and theoretically. The experimental data 
is particularly sparse for a large group (46) of binary in- 
termetallics believed to be phase separating (see Figure 
[1]). Ab initio calculations on these systems are scarce 
compared to bec and fee systems. Cluster expansion 
[22I |23| studies of hep systems 24 -[27 1 are also relatively 
rare compared to cubic systems, probably due to the 
inherent additional computational difficulties in dealing 
with a lattice having a basis of symmetrically identical 
atoms (a general cluster expansion code for arbitrary par- 
ent lattices was not available until 2002 [28|). 



METHODS 

Computations of low temperature stability of the hep 
metallic binary systems were carried out using the HT 
framework AFLOW 14 , 2|| and ab initio calculations of 
the energies with the VASP software [13] with projec- 
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tor augmented waves (PAW) pseudopotentials 3l| and 
the exchange-correlation functionals parameterized by 
Perdew, Burke and Ernzerhof for the generalized gra- 
dient approximation 32j. The energies were calculated 
at zero temperature and pressure, with spin polariza- 
tion and without zero-point motion or lattice vibrations. 
All crystal structures were fully relaxed (cell volume and 
shape and the basis atom coordinates inside the cell). 
Numerical convergence to about 1 meV/atom was en- 
sured by a high energy cutoff (30% higher than the max- 
imum cutoff of both potentials) and dense 6000 k-point 
Monkhorst-Pack meshes [33| . 

For each system, we calculated the energies of the 
reported crystal structures 34, [35| and approximately 
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226 additional structures in the AFLOW database 29 1. 
These include the 176 configurations described in 14j, 
all the symmetrically-distinct hep-, bec-, fee- based su- 
perstructures 3f| 37 1 with up to four atoms per cell, and 
the prototypes A5, A6, A7, A8, A9, All, B20, C36, D5i 9 , 
Al 2 Zr 4 , Al 3 Zr 2 , BiHf£ [13], CdTi^uPt 7 , Cu 3 Ti2, Ga 2 Hf, 
Ga 4 Ni, Ga 3 Pt^Ga 4 Ti5, HfeSc* [H, Hf 5 Sc* [l]f, Hf 5 Pb* 
0, HfPd? [3, Hf 2 TP_Jigj], Hg 2 Pt,JTl, InTh, LiB- 
MS1/2 [laEjMoTi* (TfTMoaTi* Q NbNi 8 , NiTi 2 , 
SeTl and V 4 Zns. In Ref. [13], it was shown that the 
probability of reproducing the correct ground-state, if 
well defined and not ambiguous, is r\* c ~ 96.7% ( "reli- 
ability of the method", Eq. (3)). There is no guarantee 
that the true ground states of a system will be found 
among the common experimentally observed structures 
or among small-unit-cell derivative structures. However, 
even if it is impossible to rule out the existence of an 
unexpected ground-state, this protocol (searching many 
enumerated derivative structures and exhaustively ex- 
ploring experimentally reported structures) is expected 
to give a reasonable balance between high-throughput 
speed and scientific accuracy to determine miscibility (or 
lack thereof) in these alloys. 



RESULTS AND DISCUSSION 

Of the 105 hep-hep metallic binary systems, 46 were 
previously reported as phase-separating [341I35T ]. Our cal- 
culations show that 18 of those are actually compound- 
forming at low-temperatures. Most of the predicted com- 
pounds are not hep-derived structures. These results are 
summarized in Figure [T] where all the hep metals are 
ordered according to Pettifor's chemical scale x 0, 0]- 
The symbols in each intersecting cell indicate the phase- 
separating or compound-forming nature of the binary al- 
loys. Circles denote cases in which the HT calculation 
results match the experimental data for phase separating 
(light gray) and compound-forming (dark gray) systems. 
Alloys where discrepancies are found, i.e. alloys reported 
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FIG. 1: (Color online) Compound-formation or phase- 
separation in binary systems of the hep metals, ordered by 
the Pettifor x _sca l e - Circles denote agreement of the em- 
pirical data and HT calculations on a database of ~226 
structures, for phase separating (light gray) and compound- 
forming (dark gray) systems. Squares denote HT prediction 
of stable ordered structures in systems reported as phase- 
separating (green, red or mixed for hep, non-hep or both types 
of predicted structures, respectively). 



as phase-separating but predicted to have stable ordered 
structures, are denoted by squares, green - if the pre- 
dicted structures are hep-derived, or red - if they are not 
hep-based. Squares with both colors indicate a predic- 
tion of both types of ground-states in the same system. 
No reverse discrepancies were found, i.e. no cases were 
found where systems reported to be compound forming 
were predicted to be phase separating at low tempera- 
tures. 

Most of the 46 reported phase-separating systems are 
grouped in this map into three clusters. 10 systems in 
the upper left corner, 6 in the lower right corner and 
23 in a large cluster in the center. The last metal in 
the map, thallium, has 10 phase separating systems, 7 
of which stand out of these three clusters. The ab ini- 
tio results predict 4 compound-forming systems grouped 
in the center of the upper left cluster, which is reduced 
to just 6 systems. The central cluster is also reduced to 
just 12 systems. In particular, Mg-Zr and Hf-Mg, which 
form an extended branch of this cluster are predicted to 
be compound forming. The new smaller cluster is thus 
better localized in the chemical scale space. Alternative 
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FIG. 2: (Color online) The AB structure map for hep alloys, from empirical data (upper entry in each cell) and HT calculations 
(lower entry). Circles with a single entry indicate agreement between experiment and calculation. Predicted ordered structures 
in reported phase-separating systems are highlighted in red. Discrepancies between observed and calculated structures are 
highlighted in blue. Suspected compounds with unknown prototype are denoted by a question mark. New prototypes are 
marked by * and described in Table [T] The energy difference between reported and calculated structures or between the 
reported structure (unstable in the calculation) and a two-phase tie-line is indicated in square parentheses. 



Mendeleev number assignments are sometimes suggested 
(e.g., see Ref . Hoh . One such suggestion places Mg and Be 
adjacent to each other before Sc and Y. The map in Fig. 
[T] indicates that the ordering according to the Pettifor x 
scale locates them more appropriately at the lower part 
of the scale. A few of these alternative arrangements also 
place Sc before Y, whereas our calculations indicate that 
the x scale is again more appropriate, since it groups to- 
gether the predicted compound-forming systems in the 
upper left cluster. Two isolated phase-separating sys- 
tems, Sc-Tl and Hf-Tl, are predicted to be compound- 
forming, so it now seems that TI should be more ap- 
propriately placed between Cd and Zn. The calculation 
results thus complement the empirical Pettifor map and 
group all the phase-separating systems into three small 
well-defined clusters. 

Figures |2HH show the resultant structure maps 
(XAiXb) for the binary AB, AB2 and AB3 compounds. 
The ordered structures are represented by their Struk- 
turbericht designation or prototype. A circle indicates 
agreement of the empirical data and the HT calculations. 
Cases of disagreement are indicated by rectangles. The 
upper and lower entries in each cell give the empirical 
data and the HT result, respectively. In cases where two 
different structures are reported in the literature or found 
degenerate in the calculations, both are indicated as up- 



per or lower entries. A dash (-) indicates the absence 
of a stable structure. Unidentified suspected structures 
are denoted by a question mark (?). Some of the pre- 
dicted phases (marked by an asterisk *) have structures 
for which no prototype is known and no Strukturbericht 
designation have been given. These new prototypes are 
described in Table |U 

Systems reported as phase-separating but predicted by 
our HT results to be compound forming are highlighted 
in red. In systems reported to be compound-forming, 
discrepancies between observed and calculated structures 
are highlighted in blue. In these cases, the energy dif- 
ference between the structures or between the reported 
structure (unstable in the calculation) and the two-phase 
tie-line is indicated in square parentheses. 

Most of the predicted compounds in the AB map are 
concentrated in a large cluster at the center, in the largest 
gap of the empirical map. Seven of them form a new clus- 
ter of B19 structures. Three Mg compounds of the CdTi 
prototype neighbor the reported group of Cd compounds 
with the same structure. The main feature of the empir- 
ical AB map is a large cluster of B2 compounds. The ab 
initio calculations largely preserve this cluster and the 
separation between it and the adjacent group of CdTi 
compounds. Discrepancies between reported and pre- 
dicted structures appear mostly at the margins of these 
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FIG. 3: (Color online) The AB2 structure map for hep alloys, from empirical data and HT calculations. Notations are as in 
Figure [2] § indicates that the reported prototype relaxed to the predicted one in the calculation. 



clusters. For example, the MgSc compound, reported as 
B2, is predicted to switch to the adjacent CdTi cluster 
at low temperatures. An unidentified TcZr compound 
(denoted by a question mark) is predicted to belong to 
the B2 group and an unreported MgRu compound is pre- 
dicted to belong to the CdTi cluster. Most of the energy 
differences between reported and predicted structures are 
small. Relatively large differences are found for CoZr and 
CoHf, where the reported B2 structure is predicted to be 
replaced by degenerate B33/T1I structures, and for ReTi, 
where an unreported prototype, MoTi [14j . is predicted 
to replace B2 at low temperatures. 

The structure map of the AB2 compounds is shown 
in Figure [3J In this map there are eight predicted struc- 
tures in reported phase-separating systems. Five of them 
appear in the Sc-Zr, Hf-Zr, Hf-Ti, Hf-Tl and Be-Zn sys- 
tems for which there are none in the AB map (Figure [2]) . 
The main features of the empirical map are maintained, 
which are a large cluster of C14 compounds and a smaller 
C15 group bordering it in the lower left part and a small 
CuZr2 group on the upper left. Discrepancies between 
reported and calculated structures appear at the mar- 
gins of these clusters. In addition, small groups of C49 
and C37 compounds are predicted in the upper part. 

Figure U shows the structure map for the AB3 com- 
pounds. Predicted structures appear in the Hf-Sc and 
Re-Tc systems, for which there are none in the AB 
and AB 2 maps (Figures [5] and Thus, all the sys- 



tems reported as phase-separating but predicted to be 
compound-forming have ordered structures in at least 
one of these three maps. In addition, a large cluster of 13 
structures is predicted in the non compound-forming gap 
in the center of the map. Eight of these are DO19 struc- 
tures and two have a closely related new hep-prototype 
ReaRu (and OsRes). There are no major features in 
the AB3 empirical map except a small group of DOn Y- 
compounds in the upper part. This group is reproduced 
and extended in the calculations. The reported proto- 
type CdaEr, of CdaY, is found to be degenerate with the 
DO19 structure. 

Also predicted is a cluster of Hf and Ti-compounds 
of the unreported prototype MosTi [14|. This concen- 
tration of distorted bec structures (MosTi is an 0I8 dis- 
tortion of a bec superlattice, Table QJ in six Hf and Ti 
binary systems may reflect their strong affinity to bec en- 
vironments. At zero temperature, bcc-Ti is dynamically 
unstable havin g an imaginary frequency phonon branch 
at the N point [4(1 Hl| , but it stabilizes at 882°C, where 
Ti undergoes a hep-bec transition [34j . Hence, an intro- 
duction of an alloying agent that stabilizes this phonon 
branch, can push the Ti alloy into a bec-type environ- 
ment even at low temperature. Hf also undergoes an hep- 
bec transition at a much higher temperature of 1743°C, 
[341 ] . It is known, however, that alloying elements such as 
Os, Re and Mo reduce this transition temperature and 
increase the stability domain of the bec-phase 35] . The 
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FIG. 4: (Color online) The AB3 structure map for hep alloys, from empirical data and HT calculations. Notations are as in 
Figure [3] £ In the case of, YCd3, the high throughput results find that the CdaEr structure and DO19 have the essentially same 
energy. 



calculations show that a few alloying agents may be able to drive this extension down to very low temperatures. 



TABLE I: Geometry of new prototypes (marked by "*" in Figures [2I3I4[) . Atomic positions and unit cell parameters are fully 
relaxed. 



System 


CoY* 


Sc 2 Zr* 


MoTi*[14] 


Mo 3 Ti*[14] 


Lattice 


Rhombohcdral 


Monoclinic 


Orthorhombic 


Orthorhombic 


Space Group 


R3m #166 


C2/m #12 


fmma #74 


Immm #71 


Pearson symbol 


hR4 


mS12 


o!8 


oI8 


Prim. vect. 


(hexagonal axes) 


(SG option 1) 


(SG option 2) 




a, b, c (A) 


3.896, 3.896, 20.609 


11.800, 3.271, 7.623 


4.479,3.182,9.046 


4.444, 3.173, 8.971 


a, P, 7 (dog) 


90, 90, 120 


90, 76.744, 90 


90, 90, 90 


90, 90, 90 


Wyckoff 


Col 0,0,0.3326 (6c) 


Scl 0.1914,0,-0.1370 (4i) 


Mol 0,1/4,0.3847 (4c) 


Mol 0,0,0.2440 (4i) 


positions 


Yl 0,0,0.0869 (6c) 


Sc2 0.1334,0,0.4740 (4i) 


Til 0,3/4,-0.1343 (4c) 


Mo2 0,1/2,0 (2d) 


[42, 431 




Zrl -0.4739,0,0.1931 (4i) 




Til 1/2,0,0 (2b) 


AFLOW label [29] 


"537" 


"539" 


"543" 


"541" 


System 


RcTi 2 * 


Hf 2 Tl* 


Bc 2 Zn* 


RC3 Ru* 


Lattice 


Hexagonal 


Tetragonal 


Orthorhombic 


Orthorhombic 


Space Group 


P3ml #164 


14/mmm #139 


Fmmm #69 


lmm2 #44 


Pearson symbol 


hP3 


tI6 


oF12 


oI8 


Prim. vect. 


(SG option 2) 


(SG option 2) 


(SG option 2) 




a, b, c (A) 


4.381, 4.381, 2.835 


4.422, 4.422, 7.385 


3.780, 2.0978, 10.300 


9.005, 2.757, 4.775 


a, P, 7 (dog) 


90, 90, 120 


90, 72.577, 90 


90, 90, 90 


90, 90, 90 


Wyckoff 


Rel 0,0,0 (la) 


Hfl 0,0,0.1746 (4e) 


Bel 0,0,0.17832 (8i) 


Rel 1/4,0,0 (4c) 


positions 


Til 1/3,2/3,0.299 (2d) 


T12 0,0,1/2 (2b) 


Znl 0,0,1/2 (4b) 


Rc2 0,1/2,1/6 (2b) 


[42, 431 








Rul 0,0,2/3 (2a) 


AFLOW label [29] 


"545" 


"547" 


"549" 


"551" 



G 



CONCLUSION 

Empirical structure maps employ correlations between 
simple properties of the elements (e.g. position in the pe- 
riodic table) and known crystal types to predict stable 
ordered structures of new alloys. In contrast, ab initio 
methods do not use previously available data but try to 
determine structure by computing from scratch a quan- 
tum mechanical description of the system. The structure 
maps presented above merge the results of the empirical 
Pettifor maps and HT ab initio calculations. The HT re- 
sults reveal unsuspected clusters of structures both in bi- 
nary systems known to be compound forming and in the 
empty gaps in the empirical maps, that are due to lack of 
experimental information. It is thus demonstrated that 
this integration of the two different approaches produces 
enhanced maps that should provide a more comprehen- 
sive improved foundation for rational materials design. 
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